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ABSTRACT 

The replication of the genome is a spatio-temporally 
highly organized process. Yet, its flexibility through- 
out development suggests that this process is not 
genetically regulated. However, the mechanisms 
and chromatin modifications controlling replication 
timing are still unclear. We made use of the promin- 
ent structure and defined heterochromatic land- 
scape of pericentric regions as an example of late 
replicating constitutive heterochromatin. We mani- 
pulated the major chromatin markers of these 
regions, namely histone acetylation, DNA and his- 
tone methylation, as well as chromatin condensa- 
tion and determined the effects of these altered 
chromatin states on replication timing. Here, we 
show that manipulation of DNA and histone methy- 
lation as well as acetylation levels caused large- 
scale heterochromatin decondensation. Histone 
demethylation and the concomitant decondensa- 
tion, however, did not affect replication timing. In 
contrast, immuno-FISH and time-lapse analyses 
showed that lowering DNA methylation, as well as 
increasing histone acetylation, advanced the onset 
of heterochromatin replication. While dnmt1~^~ cells 
showed increased histone acetylation at chromo- 
centers, histone hyperacetylation did not induce 
DNA demethylation. Hence, we propose that histone 
hypoacetylation is required to maintain normal het- 
erochromatin duplication dynamics. We speculate 
that a high histone acetylation level might increase 
the firing efficiency of origins and, concomitantly. 



advances the replication timing of distinct genomic 
regions. 

INTRODUCTION 

The duplication of the genome prior to ceU division occurs 
in a spatially and temporally organized manner (1). 
Studies describing DNA replication on a cytological 
level have demonstrated that sites of DNA synthesis are 
detectable as discrete nuclear foci, each corresponding to a 
cluster of several simultaneously active replication forks 
(2). On a spatial level, the number and location of these 
foci change as the cell progresses through S-phase, result- 
ing in distinct nuclear patterns associated with early, 
mid, and late rephcating regions (3-6). The temporal 
order of replication reflects the higher order organiza- 
tion of the genome. During the first half of S-phase, eu- 
chromatic regions are replicated followed by facultative 
heterochromatin during mid S-phase and finally consti- 
tutive heterochromatic regions in the second half of 
S-phase (7,8). 

The control of this spatial and temporal organization is 
intrinsically related to the coordinated firing of replication 
origins at distinct chromatin regions (9). Understanding 
the nature of this regulation has been complicated by the 
fact that a consensus sequence for rephcation origins has 
not been identified in higher eukaryotes (10), hampering 
the studying of replication origin firing dynamics. 
Moreover, certain regions change their replication timing 
during development or according to gene activity (11-13). 
This observed plasticity has led to the suggestion that rep- 
lication timing is not sequence driven, but rather modu- 
lated by dynamic changes in chromatin structure and 
composition (14). Modifications of DNA and histones 
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have been demonstrated to play a central role in defining 
chromatin structure on a local level, making them ideal 
candidates for regulators of replication timing (15). 

Several Hnes of evidence support the idea that replica- 
tion timing is regulated by chromatin structure. Specific 
chromatin modifications are known to correlate with rep- 
lication timing, such as histone acetylation with early rep- 
lication in Drosophila (16). For example, changes from 
early to late replication are accompanied by repackaging 
of DNA into nucleosonies containing deacetylated 
histones (17). In contrast, other chromatin markers, 
including DNA and H3K9 (tri)methylation (H3K9m3), 
are associated with late rephcation at particular loci 
(18). Disrupting the normal pattern of chromatin modifi- 
cations within the nucleus has been associated with 
changes in rephcation timing. For instance, inhibition of 
histone deacetylases (HDACs) by treatment with tri- 
chostatin A (TSA) leads to earher replication of the 
P-globin locus (19). In yeast, deleting the histone deacety- 
lase Rpd3 results in earher replication of late origins (20). 
Interestingly, knocking down HPl in Drosophila Kc ceUs 
results in delayed replication of euchromatic repeats and 
advanced replication of centromeric heterochromatin (21). 
Altogether, these data indicate a connection between the 
epistate of chromatin and its replication timing. However, 
the exact nature of this relationship, especially in 
mammals, remains unclear. 

To address how chromatin modifications influence rep- 
lication timing, we focused on the replication of constitu- 
tive heterochromatin (22). Constitutive heterochromatin 
exhibits a complex chromatin landscape, marked by high 
levels of DNA methylation, H3K9m3 and histone hypo- 
acetylation. These modifications help define a highly con- 
densed nature, and in mouse ceUs these regions assemble 
into higher order aggregates known as chromocenters 
(23). These structures are composed of ~10^ major satel- 
hte DNA repeats (24) and can be visuahzed by DNA 
staining with DAPI as round, highly condensed structures. 
Due to their prominent structure, chromocenters represent 
an identifiable landmark within the nucleus that can be 
easily visualized throughout the cell cycle. It is thus pos- 
sible to directly image the interactions of the replication 
machinery with the chromocenters during S-phase. 

We have undertaken a comprehensive study investi- 
gating the role of chromatin markers in defining the late 
replication timing of constitutive heterochromatin. The 
late replication timing of chromocenters was assessed by 
quantification of rephcation patterns and colocalization of 
early replication foci with chromocenters by immuno- 
FISH. We also performed live ceU imaging using fluores- 
cently tagged proteins to label the replication machinery 
and major satelhte repeats and followed chromocenter 
replication in real time. We manipulated the epistate of 
constitutive heterochromatin by altering aU three of its 
chromatin hallmarks: histone hypoacetylation, accumula- 
tion of H3K9m3 and DNA methylation. These manipula- 
tions allowed us to identify several conditions in which 
constitutive heterochromatin underwent hyperacetylation. 
This increase was associated with an earher onset of rep- 
lication. We therefore conclude that hypoacetylation, but 
not DNA methylation or accumulation of H3K9m3, plays 



a critical role in defining late-replicating regions of the 
genome. Our results put into context observations from 
different model systems, as well as from in situ and in vitro 
experiments and demonstrate that histone hypoacetyla- 
tion is an indispensable mechanism determining late rep- 
lication timing. 

MATERIALS AND METHODS 

Cells 

Wild-type mouse fibroblasts, suv39hl/2 double null, p53~^~ 
null as well as p53~^~ and dnmtl~^~ double knockout cells 
(MEF-WT, MEF-D15, MEF-P, MEF-PM, respectively) 
were grown as previously described (25,26). CeUs used 
for immunofluorescence analysis were grown on glass 
coverslips. TSA (Sigma-Aldrich, Steinheim, Germany) 
treatment was performed with a final concentration of 
20 nM over 72 h changing medium with fresh TSA every 
day. To analyze replication patterns, cells were pulse 
labeled with 100 |iM BrdU or 5]\M EdU for 30min. 
Cells were fixed for immunofluorescence in 4% formalde- 
hyde for lOmin at RT. 

Immunofluorescence 

Cells were permeabilized with 0.5% Triton X-100 for 
20min at RT fohowed by immunostaining as pubhshed 
(27). For detection of PCNA, incubation for lOmin in 
ice-cold methanol after formaldehyde fixation was neces- 
sary. Acetylated histones were detected with rabbit 
anti-H4K8Ac antibodies (1/200, Upstate, Lake Placid, 
USA) in 4% BSA for Ih at 37°C (Figure IB). 
Incorporated BrdU was recognized by mouse anti-BrdU 
antibody (5 ng/ml, Becton Dickinson, Franklin Lakes, 
USA) in conjunction with 10|ig/(il DNase for Ih at 
37°C in 1% BSA/30mM Tris HCl (pH 8.1)/0.33mM 
MgCli/l mM mercaptoethanol. Cells were then washed 
with 0.5% BSA/lniM EDTA/PBS + 0.01% Tween to 
stop DNase digestion. EdU was detected using ClickIT 
chemistry (Invitrogen, Carlsbad, CA, USA) as described 
in (28). DNA was counterstained with DAPI for lOmin at 
RT (1 |ig/ml, Sigma-Aldrich). FoUowing secondary anti- 
bodies were used: Streptavidin-Alexa 488, 1/500 
(Invitrogen), donkey-anti-mouse IgG-Cy3 and donkey- 
anti-rabbit IgG-Cy3 1/200 (The Jackson Laboratory, 
Bar Harbor, USA), goat-anti-mouse IgG-Alexa 488 
1/200 (Invitrogen). 

Immuno-FISH 

Immuno-FISH was performed as described (29) using 
mouse-anti-PCNA (1/200, Dako, Carpinteria, USA, 
Figure 4), rabbit-anti-H4K8ac (1/200, Upstate, Cat. no. 
06-760, Figure 3) antibodies and DNA probes against 
mouse repetitive sequences that were labeled with bio- 
tinylated dUTPs by nick translation. FISH signal was 
detected using Streptavidin-Cy5 1/500 (Amersham 
Biosciences, Piscataway, USA). 
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Figure 1. Manipulation of chromatin marlcs and organization of constitutive heterocliromatin. Global histone acetylation levels were assessed by (A) 
western blot analysis of H3, H3ac, H3K9ac, H4, H4K8ac and H4K12ac, as well as (B) immunofluorescence in situ using antibodies against H4K8ac. 
Confocal mid-sections of in situ stainings were imaged and the mean value of the fluorescence signal was plotted as a ratio to control cells (MEF WT 
and /)Ji^'^, respectively). For western blot analysis, the fluorescence signal of histone modifications was double normalized to the amount of histone 
signal itself (H3 and H4, respectively), as well as to control cells. While TSA treatment resulted in a clear increase of histone acetylation, dnmtl^^^ 
showed a slight increase in the level of global histone acetylation (B). On the other hand, siiv39h]/2^'^ exhibited, as expected, decreased levels of 
H3K9m3 (A). (C) Pyrosequencing analysis after bisulfite conversion demonstrated that exclusively dnmtl^^^ had abnormal levels of DNA methylation 
at major satellites, decreased from 70% to 10%. (D) Mid-confocal sections of cells stained with DAPI were used to quantify the standard deviation of 
DAPI histograms, as a measure for the homogeneity of DNA coinpaction over the nucleus (see Supplementary Figure S5). The results, presented as a 
ratio to control cells, showed the effect of the different modifications on condensation of constitutive heterochromatin. High resolution 3D-SIM iinages 
are presented to illustrate how the disruption of all three factors, histone hypoacetylation, H3K9m3 and DNA methylation, resulted in changes in the 
structural conformation of chromocenters, with TSA having the most prominent effect, as seen by a more homogeneous DAPI staining. Error bars 
represent (A) standard deviation, (B-D) 95% CL Statistical significance was tested by two-tailed /-test. **/'<0.01; ***/"< 0.001. Scale bar: 5 |rm. 



Western blot analysis and quantification 

Cells were harvested, boiled in 1 x Laenimli sample buffer 
and analyzed on western blots using the following primary 
antibodies (Figure lA): rabbit-anti-H3 (1/5000, Upstate, 
Cat #: 07-690), rabbit-anti-H4 (1/1000, Upstate, Cat #: 
07-108), rabbit-anti-H3ac (1/500, Upstate, Cat #: 06-599), 
rabbit-anti-H3K9ac (1/500, Upstate, Cat #: 06-942), 
rabbit-anti-H4K8ac, dilution 1/1000 (Upstate, Cat #: 
06-790) and rabbit-anti-H4K12ac (1/1000, provided by 
T. Jenuwein). The following secondary antibodies were 
used: anti-rabbit IgG-Alexa Fluor 647 and anti-rabbit 
IgG-Alexa Fluor 680 (1/4000, Invitrogen). Blots were 
imaged using 700 nm excitation and quantified on a 
LI-COR Odyssey Infrared Imaging System using Odyssey 
VI. 2. 15 software (Biosciences, Lincoln). Integrated pixel 
intensity was measured for each band and the respective 
background signal was subtracted. Signals were normal- 
ized to the respective loading control (histone H3 or 
histone H4) and the fold difference to the 



respective control cells was calculated using Excel 
software (Microsoft, Redmont, USA). 

Bisulfite treatment and pyrosequencing 

Genomic DNA extraction was performed using QIAmp 
DNA extraction kit (Qiagen, Valencia, USA) according to 
the manufacturer's instructions. Bisulfite conversion of 
1.5 ^g DNA per sample was performed using EpiTect 
(Qiagen). Major satelhte repeats were amplified by poly- 
merase chain reaction using the following primers; AAAA 
TGAGAAATATTTATTTG (forward) and CCATGATT 
TTCAGTTTTCTT (reverse). Three amplification reac- 
tions were pooled together and sent for pyrosequencing 
to Varionostic (Ulm, Germany). Two pyrosequencing re- 
actions, from the 3'- and 5'-end respectively, were per- 
formed per sample in order to cover the eight CpGs on 
each major satelhte unit (Figure IC and Supplementary 
Figure SI). 
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Microscopy 

3D z-stacks of fixed cells were acquired with a Leica 
TCS SP5 coiifocal laser scanning microscope (Leica 
Microsystems, Wetzlar, Germany) using a HCX PL 
APO 63x/1.40N.A. oil objective and laser lines at 
405 nm, 488 nm, 561 nm and 633 nm. To accurately com- 
pare the treated cells between different experiments, all 
images were taken using identical settings. High-resolution 
DAPI images were acquired and reconstructed with 3D 
structured illumination microscopy as described in Ref. 
(30). Quantification of S-phase patterns was performed 
by epifluorescence microscopy using an Axiovert 200 
microscope with a 63 x/ 1.4 Plan-Apochromatic oil object- 
ive (Zeiss, Jena, Germany) equipped with a Sensicam 
12-bit CCD camera system (PCO, Kelheim, Germany). 

For additional colocalization analysis of early repUca- 
tion sites and chromocenters, 3D z-stacks of cells were 
collected using a Delta Vision Olympus wide-field micro- 
scope with a 60 x/ 1.4 Plan-Apochromatic oil objective 
(Olympus Corp., Tokyo, Japan) and post-processed by 
deconvolution (Applied Precision, Issaquah, USA). For 
this, we used a radially averaged PSF recorded on the 
microscope under comparable image conditions. 

Time-lapse microscopy 

Cells were transfected with mRFP-PCNA (31) and 
MaSat-GFP (32) using Amaxa (Lonza, Cologne, 
Germany) nucleofection (solution V, program B-032) 
and thereby plated on 35 (i optical dishes (Ibidi, Munich, 
Germany). Time lapse microscopy was performed on an 
UltraView spinning disc system (PerkinElmer, UK) 
equipped with temperature, humidity and CO2 incubation 
control (PeCon, Erbach, Germany). 3D z-stacks were 
acquired using a 63 x/ 1.4 Plan-Apochromatic oil (Zeiss) 
every 30min over up to 40 hours (Figure 5, 
Supplementary Figure S4 and Supplementary Movies 
1-2). 

Image analysis and quantification 

Fluorescence intensity histogram quantifications were per- 
formed using ImageJ (http://rsb.info.nih.gov/ij/). A ROI 
was selected around each nucleus and the fluorescence in- 
tensity histogram of each nucleus, its mean value and 
standard deviation were measured. To quantify the 
levels of histone acetylation the mean values of the histo- 
grams of approx. 25 cells per condition were averaged and 
normalized to untreated cells (Figure IB). These measure- 
ments were performed in triplicates. 

To quantify the decondensation of pericentric hetero- 
chromatin the standard deviation of DAPI histograms of 
approximately 25 cells per condition were averaged and 
normalized to the control (see Supplementary Figure S5). 
The same analysis was repeated on three independent ex- 
periments and averaged (Figure ID). 

The frequency of early versus late replicating patterns 
(Figure 2 and Supplementary Figure S2) was quantified by 
counting only early and late repUcating cells (sum of both 
equals n) and calculating the percentage of early or late 
patterns, respectively. Early replicating patterns are 
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Figure 2. Effects of histone hyperacetylation, loss of H3K9 
trimetliylation and DNA methylation on late replication. Modified nu- 
cleotides (BrdU or EdU) were given to the cells for 30min before 
fixation. Detection thereof and epifluorescence microscopy allowed 
the quantification of early versus late replication patterns. (A) An ex- 
emplary field in a control cell population, with cells going through early 
(e), late (1) S-phase, as well as not replicating (non S). (B) The distri- 
bution of early versus late replication patterns. In control cells, ~40% 
of replicating cells are going through late S. While siiv39hl/2^^^ did not 
show any significant effect on this distribution, both TSA-treated cells and 
dnmtl^^^ exhibited a significant decrease in the frequency of late patterns 
down to ~15%. Statistical significance was tested using the Fisher's exact 
test. The doubling time of the control, mutant and treated cells are shown 
below the histogram. Cells grown under starvation conditions have a 
similar doubling time as TSA treated cells, even though the replication 
pattern distribution of the former is unchanged. (C) Summarizes the chro- 
inatin changes in the different cells. Black arrows indicate differences to 
the respective control cells. Scale bar: 10|xm. 
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characterized by foci distributed homogenously through- 
out the nucleus with exception of the nuclear and nucle- 
olar periphery. Late replication patterns are clearly 
recognizable due to the appearance of fewer but larger 
clusters of a couple of hundreds replication foci (2-6). 

Colocahzation of major satelhtes and early replication 
sites (Figure 4A and Supplementary Figure S3) was 
calculated using custom written software in the Priithon 
image analysis platform (http://code.google.eom/p/ 
priithon/). Images were processed using a 3D median 
filter for chromocenters and 3D Gaussian-of-Laplace 
filter for replication foci. Filtered images were thresholded 
automatically using the Otsu algorithm (33). The thresh- 
olded images were used to calculate the colocahzation per- 
centage. For this, the number of all overlapping voxels 
was divided by the total number of voxels corresponding 
to chromocenter signals. 

The total histone acetylation signal at the 
chromocenters (Figure 3) was quantified as described 
above, but after obtaining a mask for the acetylation 
channel this was multiplied by the raw acetylation 
image. The intensity of all remaining voxels was then 
summed up. 

The rate of DNA synthesis (Figure 4B) was quantified 
by measuring the integrated intensity of newly replicated 
DNA signal throughout the nucleus in ceUs pulse-labeled 
for 15min with EdU. 

S-phase length (Figure 4C) was calculated from the per- 
centage of S-phase cells in an asynchronous population 
and the total cell cycle duration. 

Colocahzation of replication foci and major satellites on 
live cell data was assessed by the H-coefficient (Figure 5 



and Supplementary Figure S4B), using the following 
formula: 
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Figure 3. Diimtl^^^ have increased levels of histone acetylation at 
chromocenters. Co-staining of histone acetylation and chromocenters 
by immuno-FISH and 3D confocal imaging (upper panel) allowed us to 
quantify the total histone acetylation signal on the chromocenters 
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of histone acetylation in these regions. Interestingly, lowering the levels 
of DNA methylation at chromocenters by knocking out dnintl'^' had 
the same effect on histone acetylation, pointing to the fact that this 
inodification is most likely to control late replication of constitutive 
heterochromatic regions. Error bars represent 95% CI. Statistical sig- 
nificance was tested by two-tailed /-test. ***/"< 0.001. Scale bar: 5|im. 
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where Ir, and Ig, is the intensity of the channels r and g in 
the pixel /' and Np is the total number of pixels. 



RESULTS 

Manipulation of the composition of constitutive 
heterochromatin 

In order to investigate the connection between chromatin 
markers and replication timing, we used drug treatment 
and genetically modified cell fines to manipulate the main 
features of constitutive fieterocfiromatin: histone hypo- 
acetylation, accumulation of H3K9m3 and DNA methy- 
lation. We treated wild type mouse embryonic fibroblasts 
(MEF-WT) with TSA to inhibit HDACs (34), thereby 
elevating histone acetylation. To modulate the levels of 
H3K9m3 on constitutive heterochromatin, we performed 
experiments in siiv39hl/2 double knock-out cells (35), 
which lack the enzymes responsible for this modification. 
Finally, we used dnmtl~^~ ceUs (25), with low levels of 
DNA methylation. Since dnmtl~^~ somatic cells do not 
prohferate normally, these experiments were performed 
in homozygous p53~^~ cells to increase their viabihty. As 
a control for the p53~^~/dnmtl~'~ double knock-out cells, 
we used p53'^' cells (25). 

Changes in the levels of chromatin markers were assessed 
by western blot analysis of cell extracts (Figure 1 A), as well 
as immunofluorescence stainings in situ (Figure IB) using 
antibodies specific to the histone modifications of interest. 
These data confirmed that TSA treatment promoted an 
increase of histone acetylation and that suv39hl/2~^~ 
cells have a decreased level of H3K9m3. Dnmtl~^~ ceUs 
lost DNA methylation at major satelhte repeats, as shown 
by pyrosequencing analysis after bisulfite conversion 
(Figure IC and Supplementary Figure SI). In the case of 
TSA treatment and suv39hl/2~^~ ceUs, the effects on the 
chromatin modifications were specific; TSA-treatment 
only affected histone acetylation and suv39hlj2~^~ cells 
exhibited only significant alterations in histone methyla- 
tion. In contrast, dnmtl~^~ cells, in addition to a drastic 
decrease in DNA methylation, also exhibited an increase 
in global histone acetylation (Figure IB and see below). 
This result is consistent with the idea that cross-talk exists 
between markers at constitutive heterochromatin, and 
that disrupting one modification may promote alterations 
in the overaU chromatin composition and state of the 
region (36). 

Moreover, we predict that disruption of these three 
chromatin modifications may also influence the condensa- 
tion levels of chromocenters (37). Figure ID illustrates 
that chromocenters in control cells exhibit the charac- 
teristic round, condensed structure described above. 
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Colocalization of chromocenters and early replication foci 
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Figure 4. Histone hyperacetylation at heterochromatic regions increase 
their replication concomitant to euchromatin. (A) Immuno-FISH was 
performed to co-stain replication structures (by eitlier PCNA or EdU) 
and chromocenters. Early patterns were selected and imaged by 3D 
confocal microscopy, as well as deconvolution microscopy. ROIs 
were defined automatically in the chromocenter channel by the Otsu 
algorithm and used to quantify the replication signal at chromocenters. 
The percentage of chromocenter material colocalizing with early repli- 
cation foci was plotted as ratio to control cells. Both, TSA-treated and 
dnmtl^^^ cells showed a significant increase of replication of hetero- 
chromatic sequences during early S-phase. (B) Inverted graph of DNA 
synthesis rate measured by quantifying the integrated intensity of 
modified nucleotides incorporated during a 15-min pulse imaged as 
3D confocal stacks (n > 90). (C) S-phase length normalized to control 
cells showing that the length of S-phase is anti-correlated to DNA 
synthesis rate. Error bars represent 95% CI (A and B) or standard 
deviation (C). Statistical significance was tested by two-tailed ^-test. 
**P<0.01; ***/>< 0.001. Scale bar: Sum. 



In contrast, this pattern is disrupted after TSA treatment, as 
well as in suv39lil/2~^~ and dnmtl~^~ cells. Chromocenters 
in these cells appeared more open and lacked the distinct, 
highly condensed appearance of control cells. Since the 
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Figure 5. Histone acetylation causes earlier onset of replication of con- 
stitutive heterochromatic regions. Progression of S-phase was followed 
by time-lapse microscopy of living cells transfected with niRFP-PCNA, 
marking sites of active replication, and MaSat-GFP, labeling major 
satellite repeats. The temporal information gained by this approach 
allowed us to unequivocally select cells going through early S-phase 
and quantify the colocalization of constitutive heterochromatic 
regions and replication structures prior to the onset of late replication. 
(A) Hyperacetylated chromosomes showed a significantly increased 
colocalization with PCNA during early S-phase when compared to un- 
treated wild-type cells. (B) Exemplary images of wild type cells going 
through early S-phase. The TSA-treated cell (bottom row) clearly 
shows colocalization of chromocenters and TSA. Error bars represent 
95% CI. Statistical significance was tested by two-tailed /-test. 
***/*< 0.001. Scale bars: 5 urn. 



decondensation of constitutive heterochromatin results 
in a more homogenous DNA staining throughout the 
nucleus, we measured the standard deviation of the re- 
spective DAPI histograms to quantify the degree of 
decondensation in treated and mutant cells (Figure ID). 
Since all three treatments influence chromatin condensa- 
tion, any differences in the effects on replication timing 
resulting from TSA treatment or knocking out dnmtl or 
siiv39h]/2 would exclude condensation as a primary deter- 
minant of replication timing. 

Manipulating chromatin epistate affects late replication 
of constitutive heterochromatin 

If changes in the composition of constitutive heterochro- 
matin promote alterations in its replication timing, we 
predict that the stereotypical late replication timing ex- 
hibited by constitutive heterochromatin would be dis- 
rupted. To test this model, we quantified the percentage 
of late replication patterns in S-phase cells. Modified 
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nucleotides, either 5-bromo-2-deoxyuridine (BrdU) or 
5-etliynyl-2'-deoxyuridiiie (EdU), were added to prolife- 
rating populations of cells and incorporated into newly 
synthesized DNA for 30min prior to fixation. Detection 
of these nucleotides allowed us to quantify the percentage 
of replicating cells going through late S-phase, as defined 
by the characteristic pattern formed by large replica- 
tion structures as compared to early S-phase patterns 
(Figure 2A). One caveat of this experiment is that several 
of the treatments we use increase the time required for the 
cells to complete the cell cycle. To ensure that the in- 
creased duplication time in mutant and treated cells does 
not result in a change of the replication pattern distribu- 
tion, we cultivated wild-type fibroblasts under decreased 
temperature and serum concentration (30°C, 5% FCS). 
While these cells progressed much slower through the 
cell cycle, similarly to TSA treated cells, we did not observe 
any difference in the distribution of replication patterns 
when compared to the control cells grown under standard 
conditions (Figure 2B). In an untreated wild-type popula- 
tion, ~40% of replicating cells exhibited staining patterns 
consistent with late replication (Figure 2B). While 
suv39hl/2~^~ did not show a significant change in this dis- 
tribution, in agreement with previous data (38), both 
TSA-treated cells and dnmtl~^~ cells exhibited a clear 
decrease (down to 15%) in the frequency of late replica- 
tion patterns. These results demonstrate that manipula- 
tion of either histone acetylation or DNA methylation 
leads to an alteration in the distribution of late replication 
patterns, suggesting a possible change in the replication 
timing of heterochromatic regions. 

Before investigating how these two chromatin modifica- 
tions are related to the regulation of replication timing, it 
was important to dissect the relationship between histone 
acetylation and DNA methylation at chromocenters (39). 
We showed that TSA had no effect on DNA methylation 
at heterochromatic sequences (Figure IC), while dnmtl~^~ 
cells exhibited increased levels of global histone acetyl- 
ation (Figure lA and B). It is therefore possible that the 
effect of reducing DNA methylation on replication timing 
is directly related to the accompanying effect on histone 
acetylation (40). However, since the increase in the global 
level of histone acetylation in dnmtl~^~ cells was not as 
pronounced as in TSA-treated cells (Figure lA and B), we 
refined our analysis to directly examine the levels of 
histone acetylation at heterochromatic regions. 

We used immuno-FISH to measure changes in histone 
acetylation levels specifically at chromocenters, simultan- 
eously detecting histone acetylation by immunostaining, 
and chromocenters by FISH. In contrast to the unspecific 
DNA staining by DAPl, FISH of major satelhte repeats 
clearly labels these regions based on their genetic compos- 
ition and independently of their condensation level. 
Subsequent collection of 3D confocal stacks allowed us 
to quantify the total acetylation signal at chromocenters 
(Figure 3). TSA treatment resulted in a clear increase of 
histone acetylation at heterochromatic regions with no 
change in DNA methylation. Strikingly, dnmtl~^~ cells 
showed the same degree of hyperacetylation at chromo- 
centers. Furthermore, treating dnmtl~'~ cells with TSA as 
done for WT MEF did not reduce the frequency of late 



replicating patterns (Supplementary Figure S2). The fact 
that TSA treatment of dnmtl~^~ cells does not result in 
additional effects on their replication timing indicates that 
chromocenters in dnmtl~^~ cells have lost normal histone 
hypoacetylation. It further suggests that the hyperacetyla- 
tion observed in dnmtl~^~ cells may be functionally 
equivalent to that in TSA-treated cells. 

Dnmtl~^~ cells provide a drug-free system in which 
histone acetylation is specifically increased at otherwise 
methylated regions, such as pericentromeric heterochro- 
matin. We have thus utilized two distinct approaches to 
promote the hyperacetylation of constitutive heterochro- 
matin. In both cases, we observe a decrease of late reph- 
cation patterns, raising the question of when these 
hyperacetylated heterochromatic regions are being reph- 
cated. Therefore, one possible model would be that hyper- 
acetylation of heterochromatin promotes earlier onset of 
replication. 

Histone hypoacetylation of constitutive heterochromatin is 
required for its late replication 

To test the hypothesis that hyperacetylated heterochro- 
matic sequences are replicated during early S-phase, we 
examined whether there was an increase of replication 
sites at heterochromatic regions during early S-phase. To 
this end, we used immuno-FISH to co-stain replication 
sites and chromocenters and directly examine their inter- 
action during early S-phase. In control cells, as expected, 
we detected very httle association between heterochroma- 
tin and sites of early DNA replication. In contrast, both 
TSA-treated and dnmtl~^~ cells exhibited increased 
colocalization between chromocenters and early replica- 
tion foci (Figure 4). Additionally, we analyzed the level 
of colocalization per chromocenter to assess whether the 
increased early replication of hyperacetylated heterochro- 
matic regions reflected single chromocenters being entirely 
replicated during early S-phase or a general increase of 
early replication distributed equally throughout all 
chromocenters. This analysis demonstrated that the shift 
to early replication is well distributed throughout chromo- 
centers, with a significant decrease in the percentage of 
chromocenters showing <10% overlap with early replica- 
tion foci from ~40% of all chromocenters in control cells 
to <20% of all chromocenters in both TSA-treated and 
dnmtl~^~ cells. However, we also found a small percent- 
age of chromocenters that are replicated mostly during 
early S-phase (Supplementary Figure S3). These results 
support the notion that hyperacetylated constitutive het- 
erochromatic regions are being rephcated in parallel to 
euchromatin, i.e. during early S-phase. 

To test if histone acetylation might result in faster rep- 
lication fork progression, we quantified the rate of DNA 
synthesis in TSA-treated and dnmtl"''" versus control 
cells. To this end, the total amount of modified nucleotides 
incorporated by replicating cells during a 15-min pulse 
was measured in confocal 3D stacks (Figure 4B). The 
outcome of this analysis demonstrated that in fact, 
DNA synthesis is slowed down in treated cells rather 
than accelerated, and this corresponds to their longer 
S-phase (Figure 4C). Importantly, these observations 
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rendered the increased colocalization of early replication 
foci and lieterocliromatic regions observed in TSA-treated 
and dnmtl~^~ cells even more significant. Interestingly, 
neither TSA treatment nor knocking out dnmtl resulted 
in checkpoint activation, as measured by Chkl phosphor- 
ylation, while our positive control of UV-C irradiated 
MEFs showed clear phosphorylation of Chkl (data not 
shown). These observations rule out the possibility that 
the described shift in replication dynamics would result 
from checkpoint activation. 

To further investigate the idea that increased histone 
acetylation correlates with a shift in the replication 
timing of heterochromatin, we monitored S-phase pro- 
gression in living cells, focusing again on the association 
between constitutive heterochromatin and the replication 
machinery. Cells were transfected with mRFP-PCNA, a 
central component of the rephcation machinery (3,31), 
and MaSat-GFP, a polydactyl zinc-finger protein which 
binds to major satelhte repeats (32) and is a hve-cell 
marker for pericentromeric heterochromatin. Series of 
time-lapse images were collected in 30-min intervals 
for up to 45 h (Supplementary Movies 1-2 and 
Figure S4A). Observation over such long periods aUowed 
us to distinguish between the different stages of S-phase in 
the same cell and to unequivocally identify cells in early 
S-phase. We analyzed the level of colocalization between 
replication foci and major satelhtes and observed a clear 
increase of colocalization between early replication foci 
and major satelhtes as a consequence of TSA treat- 
ment and in dnmtl ceUs (Figure 5 and Supplementary 
Figure S4). Based on these results, we conclude that 
increasing histone acetylation at constitutive heterochro- 
matic regions, results in an earlier onset of replication. We 
therefore propose that histone hypoacetylation is indis- 
pensable for maintenance of the late-replication timing 
of constitutive heterochromatin. 



DISCUSSION 

In this study, we used drug treatment and mutant cell fines 
to comprehensively assess which chromatin modifications 



are important for defining the late replicating nature of 
constitutive heterochromatin. By combining these ap- 
proaches with quantitative microscopy, we were able to 
directly investigate the connection between chromatin 
modifications and the spatial and temporal control of rep- 
lication timing. We show that treatment with TSA, as well 
as knocking out dnmtl, results in an eariier onset of rep- 
lication at chromocenters. In both cases, this effect was 
associated with an increase in the normal levels of histone 
acetylation at chromocenters (see below and Figure 6). In 
contrast, suv39hlj2~^~ cells with reduced levels of 
H3K9m3 and normal levels of histone acetylation ex- 
hibited no change in replication timing. These data there- 
fore suggest that the level of histone acetylation at a given 
genomic region is a major factor in determining its repli- 
cation timing. 

One important consideration when manipulating the 
chromatin state of any genomic region is the down- 
stream effects on chromatin structure and the binding of 
additional chromatin factors. This is particularly relevant 
with chromocenters, given that their DNA and histone 
modifications are directly related to their highly 
condensed nature and contributes to the recruitment of 
multiple proteins (41). In this case, the methods we used 
to perturb the molecular markers of constitutive hetero- 
chromatin (i.e. TSA treatment, dmntl and suv39hl/2 
knock-outs) all resulted in a large-scale decondensation 
of the chromocenters. However, only two out of the 
three cases promoted a shift in the replication timing of 
the constitutive heterochromatin. The fact that suv39hl/ 
2~'~ ceUs exhibit no clear change in their replication 
timing indicates that the chromatin decondensation we 
observe here is not directly related to an alteration in rep- 
lication timing. While it has been reported that an open 
chromatin condensation may facihtate early replication 
(42), we argue that a general decondensation on the 
scale that is visualized by DAPI staining in situ is not 
per se sufficient to promote early replication of hetero- 
chromatic regions. It should be stressed here that it is 
very likely that there are more subtle local differences in 
chromatin compaction and structure, which cannot be 



replication of constitutive lieterochromatin 




# replication foci # chromocenters 

Figure 6. Summary of the effect of chromatin changes in rephcation timing. Histone hyperacetylation at constitutive heterochromatic regions results 
in earlier replication onset. While DNA methylation has an indirect effect on replication timing via histone acetylation (dashed arrow), both 
H3K9m3 accumulation and large-scale decondensation of chromocenters are not sufficient to disrupt late replication. 
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detected by a mere DAPI staining and might affect repli- 
cation timing, possibly by increasing accessibility of, e.g. 
replication factors to DNA. Knocking out suv39hl/2 
genes also abolishes the binding of HPl to chromocenters 
(35), indicating that this important heterochromatin deter- 
minant does not contribute to the regulation of replication 
timing. We, therefore, conclude that neither H3K9m3, nor 
large-scale chromatin condensation, nor HPl binding, are 
directly involved in defining the late replication pattern 
observed for constitutive heterochromatin. Interestingly, 
knocking down HPl in Drosophila affects late replication 
of heterochromatin, demonstrating that replication timing 
in mammahan cells underlies a somewhat different control 
mechanism as in invertebrates. 

In both, TSA-treated and dnmtl~^~ cells, we observe a 
shift in replication timing, whereby normally late re- 
plicating constitutive heterochromatin is replicated 
during early S-phase. TSA treatment directly promotes 
the hyperacetylation of histones, specifically relating this 
chromatin modification to the regulation of repUcation 
timing. As DNA methylation is unaffected in TSA- 
treated cells, we conclude that this DNA modification 
does not block early replication. While DNA methylation 
loss is the most direct effect of knocking out dnmtl, our 
results (Figure 3), as well as previous studies (43), demon- 
strate that this decrease leads to elevated histone acetyl- 
ation levels at chromocenters. We therefore argue that in 
dnmtl~^~ cells, concomitant changes in the level of histone 
acetylation are involved in the shift in rephcation timing of 
constitutive heterochromatin, analogous to what we ob- 
serve in TSA-treated cells. Nonetheless, we cannot strictly 
exclude the possibility that TSA-mediated acetylation of 
non-histone proteins or that knocking out dnmtl might 
have other secondary effects, which could play a role in 
the described shift in replication timing. However, the fact 
that two different strategies to manipulate the levels of 
histone acetylation, a chemical and genetic one, result in 
similar effects on replication timing, strongly suggests that 
histone hyperacetylation is involved in promoting eadier 
replication onset of constitutive heterochromatin. Thus, 
we propose that histone hypoacetylation facilitates late 
replication timing independently of both histone methyla- 
tion and DNA methylation (Figure 6). Interestingly, our 
results show that there is no mechanism absolutely pre- 
venting early replication of heterochromatic regions. 
Indeed, we have reported that in normal untreated 
mouse cells, a certain small percentage of the usually 
late-replicating centromeric regions replicate during early 
S-phase (8). 

Ultimately, replication timing is defined by the timing of 
origin firing. The relative efficiency model of origin firing 
proposes that early origins fire more efficiently, while late 
origins have a low efficiency at the beginning of S-phase, 
which increases as S-phase progresses, thus assuring that 
potential gaps of unreplicated DNA are closed in a timely 
fashion (44). In the context of this model, our data suggest 
that histone acetylation directly or indirectly plays an im- 
portant role in defining the firing efficiency of origins and 
concomitantly the replication timing of distinct genomic 
regions. There are various processes leading to origin 
firing itself at which histone acetylation might regulate 



replication timing. For instance, it has been shown recently 
that origin firing dynamics in fission yeast can be a result 
of differences in the time of ORC binding at different 
regions (45). The binding of limiting ORC factors could 
be enhanced at acetylated regions with an open chromatin 
conformation. Origin licensing is another process that can 
be enhanced by histone acetylation, since HBO 1 -mediated 
histone acetylation in yeast has been shown to play an 
important role in the loading of the Mem 2-7 complex 
(46), necessary for origin licensing (47). An increased 
basal level of histone acetylation may therefore facilitate 
origin Hcensing. Alternatively, histone acetylation could 
also play a role in the actual firing process, potentially 
by increasing accessibility or binding affinity to limiting 
firing factors, such as yeast Cdc45 (20,48), which has been 
shown to increase the firing efficiency of inefficient origins 
(45). Further experiments using high-resolution micros- 
copy could give a more detailed insight into the struc- 
tural changes resulting from histone hyperacetylation. 
Moreover, protein-DNA interaction profiling under 
different chromatin acetylation conditions could elucidate 
how this histone modification affects binding of licensing/ 
firing factors to origins and how chromatin state influ- 
ences genome metabolism. 
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